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Abstract. In this report, I summarize the experimental results on bulk properties 
and flow presented at Quark Matter 2004. It is organized in four sections: 1) Initial 
condition and stopping; 2) Particle spectra and freeze-outs; 3) Anisotropic flow; 4) 
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, Outlook for future measurements. 

A wealthy set of data related to bulk properties and flow had been presented at this 
conference. I will highlight on the results related to initial condition, partonic collectivity 
■^j- . and thermalization in relativistic heavy ion collisions, emphasizing on the exciting 
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■rj- ■ in this report are to illustrate the discussions and I refer the readers to presentations by 
O 



various collaborations for complete dataset. Other important and exciting results, such 
as strangeness (K/tt) excitation function from SPS [I] and pentaquark searches from 
^ . various experiments [2 , will not be covered in this summary. 

■ 1. Initial conditions 

C 

One of the exciting results in forward rapidity (deuteron side) at RHIC is the observation 
of suppression of particle yields at and below intermediate pt in d+Au collisions. The 
observations [SI EJ El Ej are consistent with the qualitative derivation from the theory 
of Color Glass Condensate (CGC). Lively discussions had been ongoing about whether 
models with conventional stopping and/or shadowing can explain the effects IHj- 
Conclusive evidences may have to wait for more quantitative predictions from theory 
and experimental observations of more universal scalings related to the unique physical 
quantity of saturation momentum scale (Q s )- PHOBOS Collaboration has measured 
the charged particle multiplicity distributions in the full pseudo-rapidity (—5 < r\ < 5) 
range in p+p collisions and several centralities of the d+Au and Au+Au collisions 
BRAHMS, PHENIX and STAR collaborations have presented results in various rj 
coverage as well. The results appear to be consistent with HIJING/AMPT models, 
limiting fragmentation and in disagreement with an earlier prediction by saturation 
models [3 El- Updated predictions from saturation model agree better with the 
experimental results in the d+Au centrality dependence BRAHMS presented the 
identified particle yield as function of rapidity and transverse momentum • From the 
produced particle yields and net baryon rapidity loss, it was concluded that projectile 
Au nuclei lose 25 ± 1 TeV in the central Au+Au collisions. This is about 75% of the 
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beam energy. The rapidity loss is also a relevant factor in the interpretation of the 
suppression in the forward rapidity The tc rapidity distribution is Gaussian with width 
of a = 2.26 ± 0.02 and consistent with simple estimate from Landau hydrodynamics. 
The anisotropic flow v% as function of 77 measured by PHOBOS Collaboration also 
shows that it is quite Gaussian. These are very convincing evidences that the rapidity 
distribution of the particles (— ^energy) and its subsequent expansion is NOT Bjorken 
boost invariance. There was debate at the conference about the validation of the models 
with Bjorken boost invariance at the limit rapidity range (i.e. \y\ < 1). 



2. Particle spectra, freeze-outs, and nuclear modification factor 



The statistical model has been able to describe the yield ratios of the identified particles 
successfully [13]. The question is what it means since the same model can fit well to 
those in p+p and AA collisions ^Hj • To distinguish a thermal system from a system born 
with the thermal phase space where temperature is simply a Lagrange multiplier |14j. 
we need to measure the interactions between particles and observe the evolution of 
the system from canonical to grand canonical ensemble. The interaction of particles 
was measured through resonance yields according to their hadronic re-scattering cross 
section and collective flow of many identified particles. The evolution of the system 
from canonical to grand canonical ensemble was observed by the centrality dependence 
of the canonical suppression of strangeness in small systems. 
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Figure 1. Particle ratios in p+p and central Au+Au collisions at RHIC. The 
filled circles (open circles) are ratios in central Au+Au collisions measured by STAR 
(PHENIX) Collaboration. The triangles are particle ratios in p+p minbias collisions 
measured by STAR Collaboration. The red (light blue) bars are results from statistical 
model fit to 7r, K, p, A, 0, S and O (ir,K,p). 
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2.1. Particle yield ratios and chemical freeze- out 

Both SPS and RHIC have measured a rich spectrum of particle including 
7r, K,p, A, 0, S, in p+p and many centralities in AA collisions [Ej. Fig.^left portion) 
shows the variety of the stable particle ratios in p+p and central Au+Au collisions 
at y/s NN = 200 GeV. The filled (open) circles are measurements by STAR (PHENIX) 
in Au+Au collisions. The ratios measured by PHOBOS and BRAHMS [TT| ITH] are 
consistent with the measurements plotted in Fig. ^ The filled triangles depict STAR's 
results from p+p collisions. The red (light blue) bars correspond to statistical model fit 
to 7r, K, p, A, 0, S and ft (tt, K, p) in Au+Au collisions [THJ 120] • I list a few observations: 

• The statistical model [TH] fit to ir, K, p, A, 0, H and f2 results in T c /j = 160 ± 5, 
fiB = 24 ± 4, /i5 = 1.4 ± 1.6 and strangeness saturation factor 75 = 0.99 ± 0.07. 

• Strange particles are enhanced in Au+Au over p+p collisions ^3]- This is also 
evident from the strangeness enhancement factor as function of centrality J5| . 

• The fit to tt,K,p can not constrain the strangeness parameters to predict well the 
multistrangeness particle production. 

• The p/tt ratio from PHENIX is 40% lower than that from STAR due to different 
weak-decay feed-down correction. STAR measures inclusive p production with PDG 
definition |2~T) 12*2*] (P226 "the rates given include decay products from resonances 
with ct < 10 cm"). PHENIX's measurements are weak-decay feed-down corrected 
with A from measurements and E~ from model [2*3*] . 

The centrality dependence of S/(Npart/2) AA /S pp where S is either A, H or fl increases 
monotonically and shows no sign of saturation. However, the charged particle density 
which is more propotional to the volume of the system doesn't scale linearly with 
Npart from p+p to central Au+Au collisions. It was noticed that the A yield scales 
linearly with the charged particle density in Au+Au collisions ^Hj. The discussion 
at the conference was centered at the questions of whether the number of participants 
Npart is the right reference as a volume indicator and whether the system reaches grand 
canonical ensemble. The fact that the strangeness saturation factor is consistent with 
unity [121 EH 1201 suggests that the grand canonical ensemble is reached at least in the 
central collisions. 

2.2. Resonances as a probe of freeze- out dynamics 

Also shown in Fig. ^ are the ratios of resonance yields over stable particle yields. The 
resonances are 0, K*,A,A* with ratios of <p/K,K*/K,A/p, and A*/A. These ratios 
were also shown to have different centrality dependence [23 ESI I2E] ■ When compared to 
statistical model predictions from fitting results of stable particles, the ratios of K*/K 
and A* /A in central Au+Au collisions are low. At the same time, the <p/K ratio is very 
close to the prediction and independent of centrality while A ++ jp ratio is a little higher 
than but within the error of the prediction. We reached same conclusion by comparing 
the results from Au+Au to those from p+p collisions. This indicates that the peculiar 
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feature is model-independent. A natural explanation to these behaviors will be the re- 
scattering and regeneration of resonances via hadronic interactions. The system in the 
hadronic phase is effectively a pion gas. Resonances with short lifetime (cr ~ fm) will 
decay in an environment surrounding by pions. The daughters from these decays will 
re-scatter with pions and other particles. On the other hand, the hadrons will interact 
to generate the resonances through the reverse reaction channel of the decay. In the 
case of K*, the reaction channel will be K* < — > nK and the re-scattering channel is 
7T7T — > p. Since the nn cross section is large (o- n7T ~ Scr^x), it is a good assumption that 
the daughter pions from K* decay before kinetic freeze-out will not escape the system. 
Therefore, we will not be able to observe those K* decay in the hadron gas before kinetic 
freeze-out. The reaction channels of <fr ^ KK, K* «-> Ktc,A <-> pir, and A* «-> KA and 
their daughter rescatterings have very different cross seection and the yield ratios of 
4>/K,K*/K,A/p, and A*/A are consistent with the cross sections. We conclude: 

• Interactions between particles happen in the system. Temperature from the thermal 
fit is unlikely to be just a lagrange multiplier [TJj. 

• The balances of resonance re-scattering and re-generation are not maintained in 
the course of the expansion. The resonance yields reflect the time span of the 
expansion [2^ | 125 ^ 126] . 

• Interactions happen with finite cross section comparable to the cross section 
measured and implemented in transport models [2EII2Z1- 

• It is too idealized to assume infinite cross section in hydrodynamics all the way to 
the kinetic freeze-out. 

• It is too idealized to assume that the same chemical equilibrium in the resonance 
yields at the chemical freeze-out can be maintained up to the kinetic freeze-out. 
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Figure 2. Anti-proton spectra from p+p, dAu and Au+Au collisions at ^/s NN 200 
GeV. Squares and stars are results from p+p and d+Au, respectively. The triangles 
arc different ccntrality bins in Au+Au collisions. 
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2. 3. Spectra, radial flow, nuclear modification factor, and Cronin Effect 

The identified particle spectra have been used to extract information about radial flow 
and kinetic freeze-out. The main focus is on the issue of whether there is a single 
freeze-out and if not, are the particles with different hadronic scattering cross sections 
decoupling at different time? Detailed discussions have been presented [HI E21 EEl EE] 
with comparisons/fits of models of hydrodynamics or hydro-inspired blast wave to 
the measured identified particle spectra. Although some of the results indicated that 
hydrodynamics with one single freeze-out describes the spectra well, the others showed 
that sequential freeze-out of different particles is needed to fit the data. The discussions 
raised to the level of whether we are "abusing" (relying too much on, over-interpreting) 
the model comparisons However, the sequential freeze-out is more consistent with 
the rescattering of resonance discussed in previous section. Fig. El showed the p spectra 
from p+p, d+Au and Au+Au collisions by the four RHIC experiments [HOI EH E21 • The 
plot shows that from p+p to all centralities in Au+Au collisions, the p spectra exhibit 
binary scaling at pr~2 GeV/c. The centrality difference is prominent at low pt where 
the spectrum is flatter from p+p to central Au+Au collisions. BRAHMS showed that 
there is very little centrality dependence in the spectral shape ^2], therefore centrality 
dependence of radial flow is not as dramatic. However, this depends on the p? and 
centrality coverage. The feature of low yield at low pt and binary scaling at high p? 
is counter-intuitive to the radial flow effect since one would have expected the yields at 
high pt to increase from peripheral to central due to the flow push of low pt protons to 
higher p?. This also shows that the spectrum at Pt^2 GeV/c is not sensitive to the flow 
effect and probably should not be included in the simple blast wave fit or hydrodynamic 
model comparison. 

Suppression of high p? hadron production has been observed at RHIC in central 
Au+Au collisions relative to p+p collisions. This suppression has been interpreted as 
energy loss of the energetic partons traversing the produced hot and dense meduim. At 
intermediate p?, the degree of suppression depends on particle species. Surprisingly, at 
the range of 2^pr^5 GeV/c the baryons appear to scale approximately with number 
of binary collisions in Au+Au collisions relative to p+p collisions while the mesons are 
greatly suppressed. The new data presented in this conference tried to address the 
following questions: 

• whether the effect is a baryon- meson or a mass effect. This was achieved by 
comparing the nuclear modification factors of baryons (p, A, H, Q), light mesons 
(-7T, K, rj) to heavy mesons (0, K*). 

• whether the effect is a final state or initial state effect. This was achieved by 
comparing the nuclear modification factors of different particle species in Au+Au 
to d+Au collisions. 

• what hadronization scheme can explain the experimental observations in the most 
ecomonic way. 

Fig. El shows the Rcp (central/peripheral) of ir,K,p,K*,(f),A,E and Q presented 
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Figure 3. Rcp (central/peripheral) of tt, K, p, K *, (j>, A, S and f2. 

by STAR and PHENIX collaborations. Although higher statistics are necesary for 
particles, such as (f> and Q, it is quite convincing that the Rcp exhibits a baryon- meson 
effect from the comparison. One of the highlights in this conference is the successful 
explanation of t>2 scaling and this intriguing result from coalescence/recombination at 
the quark level [331 • Nuclear modification factors in d+Au (RdAu) (HOI EH EH EH] 
have been presented in the conference and publications as well [3j]. It showed that 
the proton (A) RdAu at the intermediate pt is larger than unity as also evident from 
Fig. [21 This contrasts to the binary scaling in Au+Au collisions. The difference 
between RdAu of pions and protons is about 20% while RauAu are different by a factor of 
~ 4 [3jJ EI] • It was concluded that the particle dependence of the suppression observed 
in Au+Au collisions is also a final state effect. The good candidate mechanism is the 
coalescence/recombination at the quark level [331 I3H] . 

3. Anisotropic flows 

The anisotropic flows measure the asymmetry of particle density in momentum space 
relative to the reaction plane. The geometry is initially asymmetric in non-central 
AA collisions. The transformation from geometric asymmetry to momentum-space 
asymmetry requires strong interactions at early stage where the geometry has not been 
blurred by the expansion. Anisotropic flow is one of the key achievements of the RHIC 
bulk program where good agreement between measured t>2 and ideal hydrodynamical 
calculations was reached [37]. The conference witness beautiful set of data including: 

• Higher harmonics: Vi,v 2 ,V4 [35] 135] HO] l4*Tj. 

• Different particle species: e, 7, vr, K, p, A, K*, 0, S, VI |21 EH S21 E2 E] ■ 

• Large coverage in rapidity and centrality: — 5 < r\ < 5 [35] 

• Non-flow, cumulant, and newer techniques [35] EH]. 
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We saw a set of v 2 measurements with not only 7r,K,p,A,K*,(j) but also multistrange 
particles 5 and Q. More importantly, the measured v 2 values as function of pt scale 
with the number of constituent quarks as shown in Fig. |3J This indicates that the 
number of the constituent quarks in a hadron is a relevant degree of freedom as in the 
case of nuclear modification function IHH] . In addition, hadrons with strange quarks 
behave the same as other particles. It has been argued that Q flow will be a signature of 
partonic collectivity jjnjEj. Now Q and 5 have large v 2 which comes from early stage of 
the interaction. At the hadron phase, hadrons interact according to the hadron-hadron 
cross section as observed from resonances. Keep in mind that all the Q excited states 
do not decay back to Q. Therefore, any violent interactions need to drive the Q flow 
at hadron stage will result in disappearance of Q. The only natural and reasonable 
explanation for the Q v 2 is that it is from partonic stage. 
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Figure 4. t>2 of identified particles as function of pr- Both axes are scaled by the 
number of constituent quarks. 



4. Outlook for future measurements 

(i) Forward rapidity. 

Looking forward is always a good idea! Results of forward spectra, jet and heavy 
flavor will be important to study the effect of gluon saturation. 

(ii) High statistics: (p,r],K*,Q Raa and v 2 . 

Run 4 at RHIC will enable us to have high statistic sample to study Raa and 

(iii) Charm flow as a probe of thermalization. 

PHENIX showed promising results from v 2 measurements of charm via its 
semileptonic decay j42 j. With direct open charm and its semileptonic decay, we 
would expect exciting results of radial flow and v 2 results from RHIC run 4 by 
STAR and PHENIX collaborations at QM05. 



Bulk properties and flow 



8 



(iv) Identified particle vi, v 2 , v 4. 

STAR showed that there are simple scalings between different harmonics 14U| . 
Look forward to seeing the results from identified particle vi,v& at QM05. 

(v) Identified particle spectra in the away-side jet. How can we use the jet quenching 
to study the thermalization of the system? 

One of the highlights at this conference is the experimental capability of 
reconstructing the soft particle spectra from a jet traversing a medium in Au+Au 
collisions at RHIC [HI- The spectral multiplicities and total scalar p? in both 
near-side and away-side increase with centrality (toward central) accompanying the 
disappearance of back-to-back correlation at higher p? of the associated particles. 
Moreover, it was shown that the away-side spectrum becomes softer from p+p 
toward central Au+Au collisions while the bulk spectrum becomes harder. We can 
roughly estimate the thermalized fraction of away-side jet from the (py) at central 
Au+Au collisions by (p T ) AuAujet = 80% (p T ) AuAu,buik + < 20%(p T )p P ,jet- Similar to 
what discussed in previous section, if the chemical composition (e.g. strangeness 
content) of the away-side soft spectra changes from canonical ensemble in p+p 
collisions to grand canonical ensemble in central Au+Au collisions, this will be a 
strong signature that the jet probe "melts" and becomes part of the thermalized 
system. One should look forward to this at QM05. 

(vi) Forward-backward rapidity asymmetry in d+Au. 

It is still an ongoing research on whether the Cronin Effect seen in d+Au collisions 
can be explained by final state effect as well jHE]- Further experimental results 
on particle dependence and forward-backward rapidity asymmetry will help resolve 
the issue |3E1E21EI3- 

(vii) Energy Scan. 

It will be very important /interesting to see whether some of the key phenomena 
have any energy dependence. 

5. Conclusions 

Intriguing physics in forward rapidity in d+Au collisions has been observed at RHIC. 
Thermal/statistical model describes well the stable particle yields up to Q at mid- 
rapidity in central Au+Au collisions. The radial flow derived from the identified 
particle spectra and the finite hadronic scattering cross section from the resonance 
yields indicate that the chemical freeze-out is followed by a kinetic freeze-out. It is clear 
that the measurements of nuclear modification function and anisotropic flow (v 2) of the 
identified particles support parton coalescence/recombination for particle production 
at the intermediate p?- The strong anisotropic flows of identified particles up to Q 
together with the observation of finite hadronic cross section from resonances are a 
strong indication of partonic collectivity. In the near future, the measurements will 
emphasize on thermalization and detailed evolution of the system. 
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